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Abstract. We have developed multi-dimensional constrained covariant density functional 
theories (MDC-CDFT) for finite nuclei in which the shape degrees of freedom 8), with even y, 
e.g., B20, 822, 830, 832, B40, etc., can be described simultaneously. The functional can be one of 
the following four forms: the meson exchange or point-coupling nucleon interactions combined 
with the non-linear or density-dependent couplings. For the pp channel, either the BCS approach 
or the Bogoliubov transformation is implemented. The MDC-CDFTs with the BCS approach for 
the pairing (in the following labelled as MDC-RMF models with RMF standing for “relativistic 
mean field”) have been applied to investigate multi-dimensional potential energy surfaces and 
the non-axial octupole Y32-correlations in N = 150 isotones. In this contribution we present 
briefly the formalism of MDC-RMF models and some results from these models. The potential 
energy surfaces with and without triaxial deformations are compared and it is found that the 
triaxiality plays an important role upon the second fission barriers of actinide nuclei. In the 
study of Y32-correlations in N = 150 isotones, it is found that, for 748Cf and 2°°Fm, {32 > 0.03 
and the energy is lowered by the {32 distortion by more than 300 keV; while for °% Cm and 
252No, the pocket with respect to 832 is quite shallow. 


1. Introduction 
The shape of a nucleus can be described by the parametrization of the nuclear surface or the 
nucleon density distribution [1,2]. One of the mostly used parametrizations is the multipole 
expansion with coefficients ),,’s; see Fig. 1 [3] for a schematic show of some typical nuclear 
shapes. Many shape degrees of freedom are important not only for the ground states or small 
amplitude collective motions, but also for large amplitude collective motions such as fission. The 
readers are referred to Ref. [4] for an overview. Here we simply mention that both the non-axial 
and the reflection symmetries should be broken in the study of nuclear ground state properties 
and potential energy surfaces. 

We have developed multi-dimensional constrained covariant density functional theories 
(MDC-CDFT) by breaking the reflection and the axial symmetries simultaneously [3-7]. In 
these theories, the nuclear shape is assumed to be invariant under the reversion of x and y axes, 
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Figure 1. (Color online) A schematic show of some typical nuclear shapes. From left to right, 
the 1st row: (a) Sphere, (b) Prolate spheroid, (c) Oblate spheroid, (d) Hexadecapole shape, and 
the second row: (e) Triaxial ellipsoid, (f) Reflection symmetric octupole shape, (g) Tetrahedron, 
(h) Reflection asymmetric octupole shape with very large quadrupole deformation and large 
hexadecapole deformation. Taken from Ref. [3]. 


i.e., the intrinsic symmetry group is V4 and all shape degrees of freedom 8), with even p, e.g., 
B29, B22, B30, B32, Bao, +++, are included self-consistently. The covariant density functional can 
be one of the following four forms: the meson exchange or point-coupling nucleon interactions 
combined with the non-linear or density-dependent couplings. For the pp channel, either the BCS 
approach or the Bogoliubov transformation is implemented. The MDC-CDFT with the BCS 
approach for the pairing is named as MDC-RMF and that with the Bogoliubov transformation 
as MDC-RHB [4]. In this contribution, we will present the formalism for MDC-RMF models 
and some results of actinide nuclei. 

In Section 2, we give the formalism of MDC-RMF models briefly. The results for actinide 
nuclei are shown and discussed in Section 3. A summary is given in Section 4. 


2. Formalism of MDC-RMF models 

In this section, we briefly give the formalism of MDC-RMF models. More details can be found 
in Ref. [4]. The starting point of a RMF model with the non-linear point coupling interactions 
is the following Lagrangian [8-16], 
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are the linear, non-linear, and derivative couplings and the Coulomb part, respectively. Mpg is 
the nucleon mass, ag, ay, ars, QTV, Bs, YS, Yv, Os; Ôv, Ors, and dry are coupling constants 
for different channels and e is the electric charge. ps, prs, jv, and jrv are the iso-scalar density, 
iso-vector density, iso-scalar current, and iso-vector current, respectively. 

Starting from the above Lagrangian and under several approximations, one can derive the 
Dirac equation for the nucleons, 


hypi = {a- p+ [Mpg + S(r)] + V(r)} di = eiti, (3) 


where the potentials V(r) and S(r) are calculated from the densities. 

An axially deformed harmonic oscillator (ADHO) basis is adopted for solving the Dirac 
equation [3-7, 17-19]. Note that a RMF model with reflection asymmetry has been developed in 
a two-center HO basis [20]. The ADHO basis are defined as the eigen solutions of the Schrödinger 
equation with an ADHO potential, 
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where w, and w, are the oscillator frequencies along and perpendicular to the z axis, respectively. 
These basis are also eigen functions of the z component of the angular momentum j, with 
eigen values K = mı + ms. For any basis state ®,(ro), the time reversal state is defined as 
Palro) = T®,(ro), where T = ioyK is the time reversal operator and K is the complex 
conjugation. Apparently we have Kg = —Ka and Ta = Ta with m = +1 being the parity. These 
basis form a complete set for expanding any two-component spinors. For a Dirac spinor with 
four components, 
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where the sum runs over all the possible combination of the quantum numbers œ = 
{nz Nr; M; Ms} and f% and g® are the expansion coefficients. In practical calculations, one 
should truncate the basis in an effective way. 


In our model the nucleus is assumed to be symmetric under the V4 group, that is, for all the 
potentials and densities we can do the Fourier series expansion, 


F(p, p, z z)= Fol P,% = + > faf P, Z J osne), (6) 


where fo(p,z) and fn(p,z) are real functions of p and z. The formalism for calculating the 
Fourier components of the potentials and densities can be found in Ref. [4]. 

Either the BCS approach or the Bogoliubov transformation has been implemented in our 
model to take into account the pairing effects. For the pairing force, we can use a delta force or 
a separable finite-range pairing force [21-23]. 

To obtain a potential energy surface, i.e., the energy of a nucleus as a function of deformation 
parameters, we make multi-dimensional constraint calculations which are equivalent to adding 
external potentials during the iteration. A modified linear constraint method was included in 
our MDC-RMF calculations [3-7]. 
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Figure 2. (Color online) Potential energy curve of 74°Pu with various self-consistent symmetries 
imposed. The solid black curve represents the calculated fission path with V4 symmetry imposed, 
the red dashed curve that with the axial symmetry (AS) imposed, the green dotted curve that 
with the reflection symmetry (RS) imposed, the violet dot-dashed line that with both symmetries 
(AS & RS) imposed. The empirical inner (outer) barrier height is taken from Ref. [24] and 
denoted by the grey square (circle). The energy is normalized with respect to the binding 
energy of the ground state. The parameter set used is PC-PK1. Taken from Ref. [5]. 


The total energy of a nucleus is obtained by substituting the densities into the expectation 
value of the Hamiltonian. The center of mass correction Eem. can be calculated either 
phenomenologically or microscopically, depending on the effective interactions. The intrinsic 
multipole moments are calculated from the vector densities by 


Qan = | Prov H)PVu(O), (7) 
where Y),,(Q) is the spherical harmonics. 


3. Results and discussions 

In Ref. [5], one- (1-d), two- (2-d), and three-dimensional (3-d) constraint calculations were 
made for the actinide nucleus 74°Pu using the MDC-RMF models with the parameter set PC- 
PK1 [25,26]. In Fig. 2 we show the 1-d potential energy curve from an oblate shape with (29 
about —0.2 to the fission configuration with 399 beyond 2.0 which are obtained from calculations 
with different self-consistent symmetries imposed: The axial (AS) or triaxial (TS) symmetries 
combined with reflection symmetric (RS) or asymmetric cases. The importance of the triaxial 
deformation on the inner barrier [27,28] and that of the octupole deformation on the outer 
barrier are clearly seen: The triaxial deformation reduces the inner barrier height by more than 
2 MeV and results in a better agreement with the empirical value [24]; the RA shape is favored 
beyond the fission isomer and lowers very much the outer fission barrier. Besides these features, 
it was found for the first time that the outer barrier is also considerably lowered by about 1 MeV 
when the triaxial deformation is allowed. In addition, a better reproduction of the empirical 
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Figure 3. (Color online) Potential energy curves of 7°4U in the second barrier region calculated 
from MDC-RMF models with the NL3* [29,30], NL-Z2 [31], PC-PK1 [25, 26], DD-ME2 [32], 
and DD-PC1 [33] parameter set, respectively. The reflection asymmetric shapes are allowed. 
The axial and non-axial symmetric results are denoted by black and red curves, respectively. 
Taken from Ref. [3]. 


barrier height can be seen for the outer barrier. It has been stressed that this feature can only 
be found when the axial and reflection symmetries are simultaneously broken [5]. 

In Ref. [5], it was also examined the parameter dependence of the influence of triaxiality on 
the second fission barrier and the lowering effect of the triaxiality on the second fission barrier 
was also observed when parameter sets other than PC-PK1 are used. In Fig. 3 we show the 
comparison of the potential energy curves of 784U in the second barrier region calculated from 
the MDC-RMF models with different parameter sets and Lagrangian forms, including meson- 
exchange ones NL3* [29, 30], NL-Z2 [31], DD-ME2 [32], and point-coupling ones PC-PK1 [25, 26] 
and DD-PC1 [33]. The results with and without triaxiality are both presented. From the figure 
we can see that, although the absolute values of the barriers differ a lot among these results, for 
all five parameter sets the second barriers are lowered by the triaxiality. The largest effect is 
around 1 MeV for NL-Z2 and the smallest one is around 300 keV for DD-PC1. Thus we conclude 
that the lowering effects of the triaxiality on the second barriers are parameter independent in 
the RMF models. 

The self-consistent three-dimensional constraint calculations are very time-consuming. From 
the benchmark calculations for 74°Pu, we learned many experiences about the important roles 
played by various shape degrees of freedom in different regions of the deformation space. For 
example, around the first fission barrier an actinide nucleus assumes triaxial and reflection 
symmetric shapes but around the second fission barrier both triaxial and octupole deformations 
are important [4]. These experiences are used in a systematic study of even-even actinide 
nuclei and the results were presented in Ref. [4]. Here in Fig. 4 we show the 1-d potential 
energy curves for even-even actinide nuclei from the fission isomer with G29 = 1.0 to the fission 
configuration with 629 = 1.7 calculated from the MDC-RMF model with the parameter set PC- 
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Figure 4. (Color online) Potential energy curves of even-even actinides nuclei in the second 
barrier regions calculated from the MDC-RMF model with the PC-PK1 parameter set [25, 26]. 
The reflection asymmetric shapes are allowed. The axial and non-axial symmetric results are 
denoted by solid and dotted curves, respectively. The binding energies are normalized with 
respect to the ground states. The empirical values are taken from Ref. [24] and denoted by full 
circles. Taken from Ref. [3]. 


PK1 [25, 26]. For comparison we present the results with and without the non-axial deformations. 
The reflection asymmetric shapes are allowed in both calculations. It was known that most of 
the outer barriers of the actinide nuclei are reflection asymmetric, which is consistent with 
the observed low-energy asymmetric fission fragments in this region. However, to estimate the 
fission barrier heights with an even higher accuracy, it is desirable to include also the effects of 
the non-axial deformations. From these potential energy curves it is clear that the non-axial 
deformations lower the outer barriers by around 0.5 to 1 MeV compared with the axial symmetric 
results. The empirical values of the fission barriers are also presented in the figure as full circles. 
Apparently, by including the triaxiality, the agreement between the theory and the experiment 
is improved. 

In Ref. [6] the non-axial reflection-asymmetric (32 shape in some transfermium nuclei with 
N = 150, namely 74°Cm, *48Cf, 7°°Fm, and 7°?No were investigated using the MDC-RMF 
model. In Fig. 5, we show potential energy curves for these N = 150 isotones. For 74°Cm, 
the ground state deformation (32 = 0.020. The potential energy curve is rather flat around 
the minimum. We denote the energy difference between the ground state and the point with 
232 = 0 by Edepth Which measures the energy gain with respect to the (32 distortion. For 246Cm, 
Eaepth is only 34 keV. For 248 Cf, 250Fm, and 72No, the minima locate at 832 = 0.037, 0.034, and 
0.025, respectively. The corresponding energy gain Egepth = 0.351, 0.328, and 0.104 MeV. As 
is discussed in Ref. [6], the occurrence of the non-axial octupole (32 correlations is mainly from 
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Figure 5. (Color online) The binding energy E (relative to the ground state) for N = 150 
isotones 746Cm (dashed line), 74Cf (solid line), 7°°Fm (long-dashed line), and ?°2No (dotted 
line) as a function of the non-axial octupole deformation parameter (32. Taken from Ref. [6]. 


a pair of neutron orbitals [734]9/2 (vj15/2) and [622]5/2 (vgg/2) which are close to the neutron 
Fermi surface and a pair of proton orbitals [521]3/2 (7 f7/2) and [633]7/2 (7713/2) which are close 
to the proton Fermi surface. 

It may be hard to conclude that these nuclei have static non-axial octupole deformations 
from these results because the potential energy curve is flat around the minimum and Egepth is 
small [6]. However, the present calculations at least indicate a strong Y32-correlation in these 
nuclei. Both the non-axial octupole parameter $32 and the energy gain Edepth reach maximal 
values at 748Cf in the four nuclei along the N = 150 isotonic chain. This is consistent with the 
analysis given in Refs. [34,35] and the experimental observation that in ?48Cf, the 27 state is 
the lowest among these nuclei [36]. 


4. Summary 

In this contribution we present the formalism and some applications of the multi-dimensional 
constrained relativistic mean field (MDC-RMF) models in which all shape degrees of freedom 
By, deformations with even u are allowed. The potential energy surfaces (curves) of 240Py and 
actinide nuclei with various symmetries are investigated. It is found that besides the octupole 
deformation, the triaxiality also plays an important role upon the second fission barriers. For 
most of even-even actinide nuclei, the triaxiality lowers the outer barrier by 0.5 ~ 1 MeV, 
accounting for about 10 ~ 20% of the barrier height. The non-axial reflection-asymmetric [32 
shape in some transfermium nuclei with N = 150, namely 24°Cm, 748Cf, 2°°Fm, and 7°2No are 
studied and rather strong non-axial octupole Y32 effects have been found in 748Cf and 7°°Fm 
which are both well deformed with large axial-quadrupole deformations, 62) ~ 0.3. We note 
that it is crucial to include the reflection asymmetric and non-axial shapes simultaneously for 
the study of potential energy surfaces and fission barriers of actinide nuclei and of nuclei in 
unknown mass regions such as superheavy nuclei. 


Acknowledgement 

This work has been supported by Major State Basic Research Development Program of China 
(Grant No. 2013CB834400), National Natural Science Foundation of China (Grant Nos. 
11121403, 11175252, 11120101005, 11211120152, and 11275248), Knowledge Innovation Project 
of Chinese Academy of Sciences (Grant Nos. KJCX2- EW-NO1 and KJCX2-YW-N32). The 
results described in this paper are obtained on the ScGrid of Supercomputing Center, Computer 
Network Information Center of Chinese Academy of Sciences. 


References 


1 
2 


3 


Tone 


Bohr A and Mottelson B R 1998 Nuclear Structure Vol I (World Scientific) 

Ring P and Schuck P 1980 The Nuclear Many-Body Problem (Springer-Verlag Berlin, Heidelberg, and New 
York) 

Lu B N 2012 Multi-dimensional constrained relativistic mean field theory and the potential energy surfaces 
and fission barriers of actinide nuclei Ph.D. thesis, Institute of Theoretical Physics, Chinese Academy of 
Sciences (in Chinese) 

Lu BN, Zhao J, Zhao E G and Zhou S G 2013 (Preprint 1304.2513) 

Lu BN, Zhao E G and Zhou S G 2012 Phys. Rev. C 85 011301(R)-5 (Preprint 1110.6769) 

Zhao J, Lu B N, Zhao E G and Zhou S G 2012 Phys. Rev. C 86 057304—4 (Preprint 1209.6567) 

Lu BN, Zhao J, Zhao E G and Zhou S G 2012 EPJ Web of Conferences 38 05003, Proceedings of the 
Int. Conf. on Nuclear Structure and Related Topics (NSRT12), Dubna, Russia, 2-7 July 2012 (Preprint 
1303.0621) 

Serot B D and Walecka J D 1986 Adv. Nucl. Phys. 16 1-327 

Reinhard P G 1989 Rep. Prog. Phys. 52 439-514 

Ring P 1996 Prog. Part. Nucl. Phys. 37 193-263 

Vretenar D, Afanasjev A, Lalazissis G and Ring P 2005 Phys. Rep. 409 101-259 

Meng J, Toki H, Zhou S G, Zhang S Q, Long W H and Geng L S 2006 Prog. Part. Nucl. Phys. 57 470-563 
(Preprint nucl-th/0508020) 

Paar N, Vretenar D and Colo G 2007 Rep. Prog. Phys. 70 691-793 

Niksic T, Vretenar D and Ring P 2011 Prog. Part. Nucl. Phys. 66 519-548 

Prassa V, Niksic T, Lalazissis G A and Vretenar D 2012 Phys. Rev. C 86 024317-12 

Meng J, Peng J, Zhang S Q and Zhao P W 2013 Frontiers of Physics 8 55-79 

Lu B N, Zhao E G and Zhou S G 2011 Phys. Rev. C 84 014328-10 (Preprint 1104.4638) 

Gambhir Y K, Ring P and Thimet A 1990 Ann. Phys. 198 132-179 

Ring P, Gambhir Y K and Lalazissis G A 1997 Comput. Phys. Commun. 105 77-97 

Geng L S, Meng J and Toki H 2007 Chin. Phys. Lett. 24 1865-1868 

Tian Y, Ma Z and Ring P 2009 Phys. Lett. B 676 44-50 

Tian Y, Ma Z and Ring P 2009 Phys. Rev. C 80 024313-7 

Niksic T, Ring P, Vretenar D, Tian Y and Ma Z 2010 Phys. Rev. C 81 054318-12 

Capote R, Herman M, Oblozinsky P, Young P, Goriely S, Belgya T, Ignatyuk A, Koning A, Hilaire S, Plujko 
V, Avrigeanu M, Bersillon O, Chadwick M, Fukahori T, Ge Z, Han Y, Kailas S, Kopecky J, Maslov V, 
Reffo G, Sin M, Soukhovitskii E and Talou P 2009 Nucl. Data Sheets 110 3107-3214 

Zhao P W, Li Z P, Yao J M and Meng J 2010 Phys. Rev. C 82 054319-14 

Zhao P W, Song L S, Sun B, Geissel H and Meng J 2012 Phys. Rev. C 86 064324-6 

Abusara H, Afanasjev A V and Ring P 2010 Phys. Rev. C 82 044303-11 

Afanasjev A V 2013 Proceedings of Fifth International Conference on Fission and Properties of Neutron-rich 
Nuclei, November 4-10, 2012, Sanibel Island, USA (Preprint 1303.1206) 

Lalazissis G A, Konig J and Ring P 1997 Phys. Rev. C 55 540-543 

Lalazissis G A, Karatzikos S, Fossion R, Arteaga D P, Afanasjev A V and Ring P 2009 Phys. Lett. B 671 
36-41 

Bender M, Rutz K, Reinhard P G, Maruhn J A and Greiner W 1999 Phys. Rev. C 60 034304—20 

Lalazissis G A, Niksic T, Vretenar D and Ring P 2005 Phys. Rev. C 71 024312-10 

Niksic T, Vretenar D and Ring P 2008 Phys. Rev. C 78 034318-19 

Chen Y S, Sun Y and Gao Z C 2008 Phys. Rev. C 77 061305(R)-5 

Jolos R V, Malov L A, Shirikova N Y and Sushkov A V 2011 J. Phys. G: Nucl. Part. Phys. 38 115103 

Robinson A P, Khoo T L, Ahmad I, Tandel S K, Kondev F G, Nakatsukasa T, Seweryniak D, Asai M, Back 
B B, Carpenter M P, Chowdhury P, Davids C N, Eeckhaudt S, Greene J P, Greenlees P T, Gros S, Heinz 
A, Herzberg R D, Janssens R V F, Jones G D, Lauritsen T, Lister C J, Peterson D, Qian J, Tandel U S, 
Wang X and Zhu S 2008 Phys. Rev. C 78 034308-6 


